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O , Abstract 
O 

We propose a new method to determine the unknown parameter associated to a self-consistent 
D ■ harmonic approximation. We check the vahdity of our technique in the context of the sine-Gordon 

. model. As a non trivial application we consider the scaling regime of the 2D Ising model away 

from the critical point and in the presence of a magnetic field h. We derive an expression that 

> 

Q> ' relates the approximate correlation length T — Tc and h. 

o 
a^ 
o 

(N 

o 



X 



Electronic address: iucci@fisica.unlp.edu.ar; naon@fisica.unlp.cdu.ai 



The so called "self consistent harmonic approximation" (SCHA) is a non-perturbative 
technique that has been extensively employed in Statistical Mechanics [0] [0 and Condensed 
Matter physics 0000 applications. Roughly speaking it amounts to replacing an exact 
action ^true by a trial action ^triai that makes the problem tractable. Usually ^triai is just a 
quadratic action that depends on certain unknown parameter Q that must be determined 
through some criterion such as the minimization of the free energy of the system. This 
approximation is intimately related the the "gaussian effective potential" [H in Quantum 
Field Theories (QFT's), a variational approximation to the effective potential which uses 
a gaussian wave functional depending on some mass parameter as the trial ground state. 
It also relies on a minimization principle often called "principle of minimal sensitivity" 
to determine the additional parameter. In this work we point out that in two-dimensional 
problems there is an alternative way to obtain the quantity Q. This method is based on 
Conformal Field Theory (CFT) Moreover, we shall show that our method yields im- 
proved results with respect to the predictions of standard SCHA in the sine-Gordon (SG) 
model and allows us to give a new description of the off-critical 2D Ising model (2DIM). In 



the former we exploit the existence of exact results [11] |12| to check the consistency of our 
proposal by obtaining a qualitatively good answer for the soliton mass. We then apply the 
same idea to the 2D Ising model at T ^ Tc and h ^ 0, a non-integrable model in which 
very few quantitative results are known [|T^ We use the fermionic representation of the 
2DIM. Since the standard SCHA is restricted to bosonic models, the new procedure is also 
an extension of the gaussian approximation to fermionic 2D theories. Our main result is an 
algebraic equation which allows to get the behavior of the correlation length as function of 
T — Tc and h. 

Let us stress that we are not introducing a new approximation but just a method to 
determine its parameter. As it is well-known the SCHA is a non controlled approximation, 
i.e. there is no perturbative parameter involved. It is then clear that the same criticism can 
be made to the present proposal. 



We shall begin by depicting the main features of the standard SCHA. One starts from a 
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partition function 
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where P/i is a generic integration measure and S'true is the exact action. An elementary 
manipulation leads to 
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for any trial action S'triai- Now, by means of the property 
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for / real, and taking natural logarithm in equation (Q), we obtain Feynman's inequality 
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In general, S'triai depends on some parameters, which are fixed by minimizing the right hand 
side of the last equation. 

At this point, in order to illustrate the procedure, we shall consider the well-known SG 
model, with action 

^true = 1-0^2 ^iP)^V{-p) + jd'x CO.{M] (5) 

where ^{p) is a scalar field and F{p) is usually of the form F{p) ~ p^. For simplicity, in this 
formula we have written the kinetic term in Fourier space but we kept the interaction term 
in coordinate space. 

As the trial action one proposes a quadratic one. 
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where VL is the trial parameter. In order to perform the standard minimization procedure, 
we first evaluate (S'true — S'triai)- The result is 



(S'true S'triai) trial ~ ^ 
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where 



and p is a normal-ordering parameter ||I6|. 

Now inserting (0) in equation (H), and extremizing the r.h.s. with respect to Q, we finally 
obtain 

^]2_^g-m/i(n)-/i(p)) (9) 

This gap equation allows to extract a finite answer for Q, depending on the mass parameter 
p (the difference — Ii{p) is finite). Note that the value of p is completely arbitrary, if 

one chooses it to be equal to the trial mass Q, the solution to the equation is 

= a. (10) 
The same result is obtained if instead of p = Q one takes p = \foL. 

Let us now present an alternative route to determine To this end we will exploit a 
quantitative prediction of conformal invariance for 2D systems in the scaling regime, away 
from the critical point. Starting from the so called 'c-theorem' [|l^ Cardy [l^ showed that 



the value of the conformal anomaly c, which characterizes the model at the critical point, and 
the second moment of the energy-density correlator in the scaling regime of the non-critical 
theory are related by 

where e is the energy- density operator, is its scaling dimension and t oc (T — Tc) is the 
coupling constant of the interaction term that takes the system away from criticality. The 
validity of this formula has been explicitly verified for several models [0 ||19|. For the SG 



model, the energy density operator is given by the cosine term, its conformal dimension is 
Ae = jS'^/AiT, t is the coupling constant a/ (3'^ and the associated free bosonic CFT has c = 1. 

Now we claim that Q can be determined in a completely different, not variational way, 
by enforcing the validity of the above conformal identity for the trial action. In other words. 
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we will demand that the following equation holds: 



/■ 1 

— / (fx\x\^ {cosp^ix) cos/3</?(0)),^i^i = — — (12) 

J Stt (2- 

which is to be viewed as an equation for the mass parameter Q. Of course, if one is interested 
in comparing the answer given by this formula with the SCHA result, when evaluating the 
left hand side of ([T^ one must adopt a regularizing prescription equivalent to the normal 
ordering implemented in the SCHA calculation. A careful computation leads to the following 
gap equation: 

"V^^ 32 n2 ^^^^ 
where we have defined the variable u = (3'^ /An {0 < u < 2) and p is the normal ordering 
parameter, as before. We see that, as in the standard SCHA equation (|]), one has different 
answers for different choices of p, but in this case, the results obtained for the values ^/a 
and Q are different. In any case one gets a non trivial dependence of Q on (3"^ in contrast 
with the SCHA. This is interesting if one recalls the physical meaning of mass gaps in the 
context of the SG model. Indeed, as it is well-known, Dashen, Hasslacher and Neveu (DHN) 



11| have computed by semiclassical techniques the mass spectrum for the SG model. It 



consists of a soliton (associated to the fermion of the Thirring model) with mass 

2 — n 

Msoi = v^, (14) 
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and a sequence of doublet bound states with masses 

2(2 -M 
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with = 1, 2, ... < {2 — u)/u. (From this last condition it is easy to see that in order to have 
bound states one must have m < 2/(A^ + 1). As a consequence there is no bound state 
for u > 1). More recently Zamolodchikov ||12|, by reinterpreting Bethe ansatz results, has 



given exact expressions for this spectrum. In particular for the soliton his formula coincides 
very well with ([l^ , except for u close to 2, where it predicts a divergence. For simplicity 
here we compare our results with equation (|14D- The first thing to note is that the masses 
in the SGM spectrum also depend on u as our prediction given by equation (|^). Thus, 
in this respect our proposal seems to be able to improve the standard gaussian prediction 



for the SGM, at least qualitatively. In order to perform a more specific and quantitative 
discussion let us compare equations (|1^) and ( p^ as functions of u. We set p = y/a, which 
corresponds to the prescription employed by DHN when deriving ([T^) and (p^Sf). The result 
is shown in Fig. 1 where one can observe a general qualitative analogy between both curves. 
In particular, for 0.7 < m < 1 {u = 1 corresponds to the free fermion point of the Thirring 
model and to the Luther- Emery point in the backscattering model [^) our prediction is in 



full agreement with the values of the soliton mass as computed by DHN. We want to stress 
that for M = 1 we get fi/ ^/a = a/3/32 ^ 0.30 whereas the value given by (JT^) is l/n ^ 0.31 
(standard SCHA yields, of course, fl/^/a = 1). 

Having checked the admissibility of our proposal in a model where exact results are 
known, it is now desirable to explore a non-trivial problem. Let us consider the 2D Ising 
model away from criticality (T ^ and /i 7^ 0): 

^ = 5*0 + j d^x[tt{x) +ha{x)], (16) 

where Sq is the critical action, t oc {T — Tc), and e(x) and cr(x) are the energy-density 
and spin operators, respectively. We shall use the fermionic representation for the above 
action. Thus Sq is a free massless Majorana action and e oc ^\E'. On the other hand, the 
expression of cr{x) in terms of the Majorana fields is more involved. Indeed, by means of a 
Jordan- Wigner transformation it can be written as an exponential of a fermionic bilinear. 
In analogy to the usual SCHA method, we propose the following quadratic trial action: 

Striai = So + n j (fxe{x), (17) 

The conformal equation (0) for the present case takes the form 

Jd'rr' [t' (2 - (e(r).(0)),„,, + h^2 - {a{r)a{0)\^,^, + 

+2th{2- A,) (2 - A.) Hr)am,,J = ^, (18) 

where we have set c = 1/2, which is the central charge corresponding to Majorana free 
fermions and A^ = 1 and A^- = 1/8 are the scaling dimensions of the corresponding opera- 
tors. Now we have to evaluate the v.e.v.'s in the trial theory. This will give us an equation 
for Q as function of t and h. The energy-energy and the energy-spin correlation functions 



have been computed by Hecht |24] whereas the spin-spin correlator can be found in the work 



of Wu, McCoy, Tracy and Barouch ||2^. As usual, one defines a correlation length ^ = 1/40 



and considers the scaling limit given by ,^ ^ oo, r oo, with r/^ fixed. The next step 
is to use the expressions of the correlators for (r/^) << 1 and perform the corresponding 
integrals. At this point we have to take into account that the correlation functions are pro- 
portional to certain scaling functions -F±(r/^) where the + and — signs correspond to the 
cases O > and O < respectively. In other words, the parameter Q can be seen as defining 
a new "effective" critical temperature, and the functions F± describe the scaling regime 
above and below this temperature. Since we are approximating a magnetic perturbation 
of the system it is clear that we must use the functions F_. Thus we obtain the following 
equation relating ^, h and t: 

t2(40' + Ci/i2(40'^/' + C2t|/i|(40'^/^ = 1 (19) 

where we have introduced the numerical constants Ci = 0.749661 and C2 = 0.186966. The 
absolute value of the magnetic field in the second term comes from the fact that (ecr) oc (cr) 
and the product (a) h has to be positive since the magnetization and the magnetic field have 
the same orientation. For ^ fixed this equation gives a simple dependence of h as function 
of t. Indeed, for h > one has a slightly rotated semi ellipse in the upper h — t plane, and 
for h < one has its reflection over the t = axis. 

An alternative form of equation ([19D is obtained if one introduces the dimensionless 

I-— 

combination X = I ^ I ^^/4^o- 
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The + and — signs in the third term of the left hand side correspond to the cases t > and 
t < 0, respectively. The action (0) defines a one parameter family of field theories which 
can be labelled by x- Moreover, the particle content of the model is expected to undergo 
drastic chang function of x [0- 



In order to check the consistency of the above equations we first consider the limits h —>■ 

and t — >■ separately. The first case corresponds to x 00 and one immediately obtains 

I— — 

^ = ^0, as expected. In the second case one has x and then we get C, ^ \ h \ , 



which is in agreement with the exact result obtained in Ref.[|2T]]. Let us mention that in 
this reference the constant of proportionality was exactly determined to be 4.4, whereas 
our approximate computation yields 3.7. Going back to the general case, we have solved 



equation (0) numerically for ^ as function of x fo^' both t > and t < 0. The results are 
plotted in Fig. 2. In the t > case the correlation length increases in a monotonous way 
from the zero value and it approaches the h = value from below as x ~^ oo. In the 
t < case, although the behaviour of ^ seems very similar to the previous case, it presents a 
subtle difference shown in Fig. 3. For x ~ 2 the correlation function goes over the value, 
reaches a maximum and then tends to from above as x ^ oo. As this behaviour depends 
on the value of the constants Ci and C2 we do not know whether this is indeed a property 
of the Ising model or an artifact introduced by our approximation. 



To conclude, we have reconsidered the well-known SCHA method in which a compar- 
atively complex action is replaced by a simpler quadratic system depending on a mass 
parameter Q which is usually determined through a variational calculation. Taking into 
account the (l+l)-dimensional case, we have proposed an alternative way for evaluating Q. 



Our proposal is based on a consequence of Zamolodchikov's [0 c-theorem first derived by 
Cardy [0]. We have illustrated the idea by considering the SG model. We showed that for 
this model our method gives a quite good prediction for the behavior of the soliton mass as 
function of (see equations (|^) and (|I4D and Fig. 1). 

As a non-trivial application we have considered the 2D Ising model away from criticality 
(T 7^ Tc and h ^ 0). Starting from a continuum field theoretical description in terms of 
Majorana fermions, we proposed a quadratic trial action depending on a parameter Q that 
defines an approximate correlation length ^. Our main result is given by equation (|19|)(or 
its alternative form (pO])) which allows to determine the parameter Q (i.e. ^) in terms of the 
original physical parameters t and h. 

It would be interesting to test our approach in other models such as the continuum version 
of the tricritical Ising model, which is described by the second model of the unitary minimal 
series [^] with central charge c = 7/10. 
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Figure caption 

Figure 1: Masses in units of ^/ol as functions of u. Tlie filled line is Mi^oij \foi-, whereas 
the dashed hne represents fi/^a as given by equation ITB. 



Figure 2: Correlation length in units of as a function of x fo^' both t > and t < 0. 



Figure 3: As Fig. 2, showing details of the behaviour of ^(x) for t > and t < 0. 



11 



